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If the unidentified emission line at ~3.55 keV previously found in spectra of nearby galaxies and galaxy clus¬ 
ters is due to radiatively decaying dark matter, one should detect the signal of comparable strength from many 
cosmic objects of different nature. By studying existing dark matter distributions in galaxy clusters we identified 
top-19 of them observed by XMM-Newton X-ray cosmic mission, and analyzed the data for the presence of the 
new line. In 8 of them, we identified > 2<r positive line-like residuals with average position 3.52±0.08 keV 
in the emitter’s frame. Their observed properties are unlikely to be explained by statistical fluctuations or as- 
trophysical emission lines; observed line position in M31 and Galactic Center makes an additional argument 
against general-type systematics. Being interpreted as decaying dark matter line, the new detections correspond 
to radiative decay lifetime r DM » (3.5 — 6) x 10 27 s consistent with previous detections. 


The origin of missing mass in cosmic objects ranging from 
dwarf galaxies to galaxy clusters, large-scale structure and the 
observable part of our Universe, remains unknown. Assuming 
Newtonian/Einsteinian gravity and dynamics to be valid one 
has to introduce new type of matter - the dark matter - pre¬ 
sumably in form of new elementary particles beyond the Stan¬ 
dard Model of particle physics. Non-zero interaction strength 
of dark matter particles with Standard Model particles may 
lead to 2-body radiative decay of the dark matter. Because 
present-day velocities of dark matter particles in haloes should 
be highly non-relativistic, such process would produce the 
narrow dark matter decay line. This motivates extensive on¬ 
going studies of such lines in spectra of cosmic objects with 
established dark matter contribution. 


public observations of European Photon Imaging Camera 
(EPIC) S 0 on-board XMM-Newton X-ray observa¬ 
tory |IS] - the most sensitive existing instrument to search 
narrow faint X-ray lines [ 16}]. 

Object selection. For distant objects, the dark matter decay 
flux [in photons cm -2 s _1 ] is 

jp _ ^DM^fov 

47T JTTdm Aim 

where <S DM = J p DM (l)dl is the dark matter column density 
along the line of sight, f!f ov <C 1 is the Field-of-View of the 
instrument, to dm - mass of the dark matter particle, r DM - 
radiative dark matter decay lifetime. 


The most intriguing dark matter decay line candidate reported 
so far is the ~3.55 keV line detected in central part of the 
Perseus galaxy cluster and different combinations of galaxy 
clusters 111, Andromeda galaxy and Perseus galaxy cluster 
outskirts y]], Galactic Center yfl, Perseus, Coma, and Ophi- 
uchus galaxy clusters a 1 , see also a. Although decaying 
dark matter can naturally explain basic properties of the line 
in these objects - correct scaling of the line signal with object 
redshift and projected dark matter mass density - several alter¬ 
native hypotheses invoking standard or anomalously enhanced 
astrophysical line emission IfUU], decay of excited dark mat¬ 
ter states & annihilating dark matter lIToL [Till , dark matter 
decaying into axion-like particles with further conversion to 
photons in magnetic field |[l2] have been proposed thereafter, 
see H and references therein. 


We define signal-to-noise ratio as 


SNR = 


Nline 

AN ba ck ’ 


( 2 ) 


where Nu ne ex S DM t 0 bs is the number of counts expected 
from decaying dark matter during observation time t 0 b s and 
A Nback is the uncertainty of background counts. Neglecting 
systematical errors and assuming Gaussian statistics, we ob¬ 
tained ANback = VBlobs, where B is background count rate 
(in cts/s) measured in 3.4-3.65 keV in object’s rest frame. So 
the signal-to-noise ratio is proportional to 


SNR oc S DM 



(3) 


To further check the decaying dark matter origin of the 
~3.55 keV line, we identified cosmic targets having the 
largest expected decaying dark matter signal. We used 


1 Based on observed ~ 60 — 100 eV shift between line positions, fl found 

improbable that their detections in Coma and Ophiuchus are of the same 

origin as in Perseus. 


Because of strong dependence of SNR from S DM , we first 
identified objects with the largest column density in XMM- 
Newton/EPlC Field-of-View. In this paper, we concentrated 
on galaxy clusters - the objects having the largest col¬ 
umn densities inside the characheristic radius of dark mat¬ 
ter haloes mill. For each of these objects, we calcu¬ 
lated dark matter column densities S 0 bj inside the innermost 
14’ radius circle Ii\.\ (roughly corresponding to the XMM- 
NewtonfEPlC Field-of-View) using dark matter distributions 
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available in & see Appendix[A] and broadened the obtained 
column density ranges by 0.15 dex to account typical residual 
uncertainties in dark matter distributions 0. We then calcu¬ 
lated the foreground column density Smw from Milky Way 
halo by using the newest dark matter distributions compiled 
in 0 . Because in most of objects XMM-Newton/EPIC energy 
resolution is comparable to the energy split between the ex¬ 
pected dark matter decay signals from object and the Milky 
Way halo, we multiplied the latter by a correction factor: 


5dm — S 0 bj + Smw x exp I — ■ 


2 rp2 
obj^line 

2fJ instr 


(4) 


where z 0 bj is the object redshift, Eu ne ss 3.55 keV is the line 
position, cTinstr is the Gaussian dispersion width correspond¬ 
ing to energy resolution of the instrument 2 . According to 
Fig. 5.24 of lUff], for XMM-Newton/EPIC imaging spectrom¬ 
eters Uinstr ~ 60 eV at the energies of our interest. 

We identified 20 galaxy clusters with the largest ,S' DM , one of 
them - Abell 539 - was not observed by XMM-Newton. The 
basic properties of the remaining 19 objects are summarized 
in Table ID 


Data reduction. For objects of our interest, we first down¬ 
loaded all public observation data files for MOS 0. and 
PN 114,] cameras of XMM-Newton X-ray observatory | Ff5il . 
and processed them using Extended Sources Analysis Soft¬ 
ware (ESAS) package [20] publicly available as part of XMM- 
Newton Science Analysis System (SAS) v. 14.0.0. Time in¬ 
tervals affected by highly variable background component 
- soft proton flares [20] - are filtered using ESAS scripts 
mos-filter and pn-filter. We used the standard fil¬ 
ters and cuts provided by ESAS software. We excluded 
bright point sources detected with the standard SAS proce¬ 
dure edetect_chain, extracted source spectra and pro¬ 
duced response matrices inside the 14’ radius circle around the 
NASA Extragalactic Database (NED) source position using 
ESAS procedures mo s- spectra and pn-spectra. Back¬ 
ground spectra were prepared by ESAS scripts mos_back 
and pnJoack. For PN camera, we additionally corrected 
the obtained spectra for out-of-time events. Finally, for each 
object we combined spectra and response files from MOS 
and PN cameras using addspec FTOOL procedure similar 
to Hi, and grouped the obtained spectra by 60 eV per en¬ 
ergy bin to make the bins roughly statistically independent. 

Spectral modeling. For each object, we modeled separately 
its combined MOS and PN spectra in Xspec spectral pack¬ 
age with the sum of non-thermal (powerlaw) and ther¬ 
mal (line-free apec) continuum components, and several 
narrow zgaussian lines of astrophysical origin and the 
new line absorbed with phabs model and folded with re¬ 


2 For Gaussian line, the energy resolution is characterized by full width at 
half-maximum (FWHM) of the line which is 2<r f 2 log(2) Ri 2.35<7j ns t r . 


sponse files. The powerlaw index F = 1.41 and nor¬ 
malization 11.6 photons cm -2 s -1 sr _1 keV -1 at 1 keV are 
fixed to best-fit values of 0. We chose the modeled en¬ 
ergy range 2.1-6.0 keV avoiding strong emission lines. To 
account residual line-like calibration uncertainties we added 
0.4% (MOS) and 0.25% (PN) systematic errors in quadra¬ 
tures using Xspec parameter systematic, according to 
Sec. 5.3.5 of 111911 . The absorption hydrogen column den¬ 
sity was fixed at weighted Leiden-Argentine-Bonn survey I 22 I 
value obtained through n H tool of the NASA High Energy 
Astrophysics Science Archive Research Center (HEASARC). 
The redshifts of apec and zgaussians were fixed at val¬ 
ues from NASA Extragalactic Database (NED). The new line 
position is allowed to vary in 3.35-3.70 keV. To calculate the 
contribution of other astrophysical lines in this region, we 
used the bright ‘reference’ S XV line complex at ~2.63 keV 
detected in majority of our combined spectra, see AppendixlBl 
for details. 

Fit quality, plasma temperatures, maximal expected fluxes 
from K XVIII line complex at 3.51 keV, new line positions 
and normalizations, and increase of \ 2 statistics due to the 
new line are summarized in Table [HI If no line is detected at 
Ict level, we put 2 ct upper bounds instead. 

Discussion. Assuming decaying dark matter origin of the line 
at ~3.55 keV previously reported by M, we identified 19 
galaxy clusters with the largest expected significance of dark 
matter decay signal. Using publicly available XMM-Newton 
observations of their central parts, we confirmed previous de¬ 
tections in Perseus Dli and Coma 0] clusters, and found 
>2 ct positive line-like residuals in 6 new objects, see Table HU 
for details. We consider the following traditional origins of 
new line detections: (a) pure statistical fluctuations; (b) con¬ 
tribution from nearby astrophysical emission lines; (c) (un¬ 
known) systematical effect. 

To check whether pure statistical fluctuations may cause 
our detections, we simulated cluster spectra using FTOOL 
fake it for all objects shown in Tablellllbased on our best-fit 
models without adding a new line in 3.40-3.65 keV, looked for 
A\ 2 increase caused by adding the new narrow zgaussian 
line in the energy range 3.40-3.65 keV (thus accounting for 
the look-elsewhere effect). The average value of 3 maximal 
A\ 2 for each simulation is in the range 2.1-5.2, much smaller 
than 13.6 (12.5) obtained from our MOS (PN) observations. 
Therefore, we conclude that pure statistical fluctuations alone 
can not be responsible for line detections in Table UTI 

Explanation of the new lines with astrophysical line contribu¬ 
tion is also unlikely. The maximal contribution of the most 
promising astrophysical line candidate - K XVIII line com¬ 
plex at ~3.51 keV - is already included to our model, see Ta¬ 
ble [II] Other astrophysical lines are both too faint and should 
produce detectable signatures at other energies. For exam¬ 
ple, to explain the excess in Virgo cluster (also consistent with 
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TABLE I: Galaxy clusters observed by XMM-Newton ranged by expected significance of decaying dark matter signal from their central parts. 


Object 

(1) 

redshift 

(2) 

S D „, M/pc 2 
(3) 

XMM-Newton ObsID 

(4) 

MOS/PN exposure, ks 

(5) 

MOS/PN SNR , arb. units 

(6) 

Virgo 

0.0036 

624-1338 

0108260201, 

0110930701, 0210270101 

193.2/62.9 

3.1-6.7 / 2.3-5.0 

Centaurus 

0.0114 

818-1721 

0046340101, 

0406200101 

348.2/ 123.7 

2.9-6.0 / 2.2-4.6 

Abell 85 

0.0551 

130-777 

0065140101, 

0723802101, 0723802201 

401.6/ 136.7 

0.6-3.8 / 0.5-2.7 

Abell 478 

0.0881 

83-969 

0109880101 


111.8/43.2 

0.3-3.5/0.2-2.7 

Abell 2199 

0.0302 

123-826 

0008030201, 

0008030301, 0008030601, 

254.4/97.5 

0.5-3.4 / 0.4-2.6 




0723801101, 

0723801201 



Abell 496 

0.0329 

124-772 

0135120201, 

0506260301, 0506260401 

250.6/81.0 

0.5-3.1 /0.4-2.2 

2A0335+096 

0.0363 

75-608 

0109870101, 

0109870201, 0147800201 

225.2 / 89.0 

0.4-2.9 / 0.3-2.4 

Abell 1060 

0.0126 

451-1420 

0206230101 


66.7 / 24.7 

0.9-2.8 / 0.8-2.4 

Abell 3266 

0.0589 

385-768 

0105260701, 

0105260801, 0105260901, 

179.8/63.9 

1.4-2.8 / 1.1-2.2 




0105261001, 

0105261101, 0105262101, 






0105262201, 

0105262501 



Abell S805 

0.0139 

286-660 

0405550401, 

0694610101 

92.2/12.5 

1.2-2.7/0.5-1.1 

Coma 

0.0231 

191-1193 

0124711401, 

0153750101, 0300530101, 

343.8 / 122.0 

0.4-2.6/0.3-2.1 




0300530301, 

0300530401, 0300530501, 






0300530601, 

0300530701 



Abell S239 

0.0635 

256-553 

0501110201 


81.0/28.1 

0.9-1.9 / 0.6-1.3 

Abell 2142 

0.0909 

88-573 

0674560201 


104.8/40.9 

0.3-1.7/0.2-1.3 

Abell 2319 

0.0557 

359-716 

0302150101, 

0302150201, 0600040101 

159.7/60.7 

0.8-1.5 / 0.6-1.2 

Abell 1795 

0.0625 

83-589 

0097820101 


71.7/23.2 

0.2-1.5 / 0.2-1.1 

Abell 209 

0.2060 

67-500 

0084230301 


33.5/11.2 

0.2-1.4 / 0.1-1.0 

Perseus 

0.0179 

418-871 

0085110101, 

0305780101 

316.4/44.2 

0.7-1.4 / 0.3-0.7 

PKS0745-191 

0.1028 

59-458 

0105870101 


31.9/5.2 

0.1-0.8 / 0.1-0.5 

Triangulum 

0.0510 

379-757 

0093620101 


18.6/ — 

0.3-0.6 / — 



3.0 3.2 3.4 3.6 3.8 4.0 

Energy [keV] 



3.0 3.2 3.4 3.6 3.8 4.0 

Energy [keV] 


FIG. 1: Examples of spectral dataset with identified extra line, see Table ITT1 for details. The spectra are binned by 60 eV and presented in 
detector’s frame similar to Qjj]. Blue and red residuals (bottom) are shown with respect to the best-fit model with and without adding an extra 
line, respectively. Left: MOS spectrum of Abell 2199. Right: PN spectrum of Abell 496. 


pure statistical fluctuation), one should assume 3.398 keV as- 
trophysical line from S XVI ^5 times higher than the maximal 
contribution from the bright S XV line complex at ~2.63 keV 
obtained from Fig. [4] On the other hand, one cannot ex¬ 


clude the possibility of strongly super-solar abundance be¬ 
cause there can be variations of Potassium abudance up to 1 
dex iHQ- According to l26h . further studies of the new line 
using forthcoming observations of Soft X-ray spectrometer 
on-board Astro-H X-ray observatory 0 of Micro-X sound- 
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TABLE II: Model parameters of MOS/PN combined spectra of galaxy clusters listed in the previous Table. Line positions are given in cluster’s 
rest frame. Column (4) shows our estimate on maximal K XVIII line flux at 3.51 keV using prominent S XVI line complex at 2.63 keV, see 
Appendix [B] for details. Errors on line position and flux are given at la level for 2 d.o.f. calculated using A\ 2 = 2.3. Line fluxes are in 
10~ 6 ph cm -2 s“\ abundances are in Solar values given by [23]]. The new line is detected at > 2 a (corresponding to Ax 2 > 6.2) in 8 
objects (marked in bold), confirming previous detections in Perseus m and Coma 0]. 


No Object 

X 2 /d.o.f. 

T e ,unei keV 

max K flux at 3.51 keV 

New line position, keV 

New line flux 

^XZine 


(1) 

(2) 

(3) 

(4) 


(5) 

(6) 

(7) 

1 

Virgo 

68.9/48 / 60.0/43 

1.4/ 1.4 

< 0.9 / < 

0.7 

3.38++/- 

4.2t^ / < 9.3 

3.8/0.4 

2 

Centaurus 

64.1/47 / 64.3/46 

2.2/2.2 

< 5.8 / < 

5.6 

3.5lt°;i 2 /- 

25. 2 + 24 ^ / < 15.6 

2.9/0.1 

3 

Abell 85 

37.5/49/61.9/46 

2.2/3.4 

<0.7 /< 

0.9 

3 44+ 0 - 06 /_ 

^•^-0.05 ' 

6.3+|^ / < 4.2 

7.0/0.0 

4 

Abell 478 

54.1/47/48.9/49 

2.2/ 1.4 

< 0.4/< 

0.5 

—/— 

< 20.4/< 13.6 

0.1/0.4 

5 

Abell 2199 

46.8/47 / 70.6/52 

2.7/2.7 

< 1.8 / < 

1.6 

3 41+ 0 - 04 /_ 

-’•‘♦A—0.04 ' 

Idling / < 10.0 

10 . 2 / 0.1 

6 

Abell 496 

36.3/45 / 66.7/48 

3.4/2.2 

< 1.3/< 

0.5 

3 ffff+0.06 / 4 C+O.O 4 

j, ‘ 5 *’-0.09 ' ^*^-0.03 

7.5i®;J /16.8+4 

6.2/18.8 

7 

2A0335+096 

70.9/49 / 65.3/49 

2.2/2.7 

< 1.5/< 

1.3 

—/ — 

< 15.5/< 10.7 

0.5/1.6 

8 

Abell 1060 

67.4/48 / 63.0/50 

2.2/2.7 

< 2.9 / < 

1.8 

-/- 

< 27.1/< 21.2 

0 . 2 / 0.0 

9 

Abell 3266 

40.7/47 / 67.2/50 

1.7/ 1.7 

< 0.3 / < 

0.3 

3 64+ 0 - 05 / 3 53 + 0.04 
-i-u^-O.OS ' ^*^-0.06 

6.5++8.71+ 

3.9/8.0 

10 

Abell S805 

49.0/45 / 33.4/26 

1.7/ 1.4 

< 0.2 /< 

0.3 

_/ 3 /:3+0.05 

/ J.OJ_ 0> 06 

< 8.7/17.11+ 

0.3/10.8 

11 

Coma 

41.2/37/54.7/48 

2.2/4.3 

< 1.9 / < 

2.0 

3 40 + 0.04 / 0 41 + 0.11 
*^•^-0.05 ' ~’*^ i - 0.10 

23.7+V / 14.81+ 

16.6/3.5 

12 

Abell S239 

56.7/48 / 60.8/52 

1.4/ 1.7 

< 0.1 / < 

0.2 

—/ — 

< 12.3/< 13.6 

0.3/0.5 

13 

Abell 2142 

63.9/50 / 56.9/50 

1.4/ 1.4 

< 0.3/< 

0.3 

— / — 

< 9.8 / < 17.4 

0 . 0 / 0.8 

14 

Abell 2319 

49.4/47/61.6/51 

1.4/2.2 

< 0.4/< 

1.4 

3 5Q+0-05 / 0 ci+0.11 
J.*?^_0.06 ' 

i8.6l 7 ° 4 7 / 10 . 5 I i 2 :® 

13.9/2.4 

15 

Abell 1795 

61.5/51/64.6/50 

1.7/ 1.7 

< 0.3/< 

0.5 

—/ — 

< 12.4/< 16.5 

0.7/0.0 

16 

Abell 209 

62.3/50 / 68.0/48 

1.4/ 1.4 

< 0.5 / < 

0.2 

— / — 

< 17.4/< 9.4 

0 .6/0.0 

17 

Perseus 

69.6/48 / 81.2/47 

2.7/2.7 

< 4.5 / < 

6.1 

3 ce+0.05 / 
J.3O-0.08 ' 

25.2lj|| / < 70.4 

9.8/0.7 

18 

PKS0745-191 

68.9/47 / 56.0/53 

2.2/1.4 

< 0.9/< 

1.5 

3 63 +0 07 /_ 

J.OJ—o 23 ' 

12.5++ / < 40.7 

2.4/ 1.6 

19 

Triangulum 

56.7/49/ — 

2 .2/ — 

< 1.4/- 


—/ — 

< 47.1/ — 

0.7/ — 


ing rocket experiment [28|] with superior spectral resolution 
< 4 — 7 eV can reveal its astrophysical origin. 

The systematic origin of the new line is shown unlikely in 
pioneering papers UJ, |2|]. In addition, we plotted in Fig. [2] 
the dependence of the line position from the object’s red- 
shift. If the new line were due to systematic effects, one 
would expect the corresponding new line in nearby (z = 0 ) 
objects at ~3.40 keV, in apparent tension with observations 
(3.53±0.03keV for M31 and 3.539±0.011 keV for Milky 
Way H) which in turn are better consistent with the new line 
generation in cosmic objects. The mean value of the line po¬ 
sitions in Fig.[2]is 3.52 keV. The average spread between line 
positions is 75 eV close to cri nstr ~ 60 eV and consistent with 
our simulations, according to that the position of ~ 3(7 line 
can be recovered with ±110 eV precision in 90% of cases. 

Interpreting the new line due to decaying dark matter 0 gives 
the radiative decay lifetime r DM « (3 — 6 ) x 10 2 ' s consistent 
with previous detections [ l-4j |29tl . see Fig. Q] Non-detection 
of the line in some of our galaxy clusters does not exclude the 
dark matter line origin; the strongest 2 cr upper bound for our 
objects comes from Virgo cluster: t dm > 3.5 x 10 27 s. 

Non-detection of ~3.55 keV line in stacked dSphs by (55] 
is also mildly consistent with these results; planned observa- 



NED redshift 

FIG. 2: Position of new line detections in cluster’s rest frame as 
function of redshift. Only detections at > 2a (corresponding to 
A \ 2 > 6.2 for 2 d.o.f.) are shown. Red and black dashed lines 
show expected behavior in case of purely systematic and cosmic line 
origins (assuming line position 3.52 keV in detector’s frame expected 
from , respectively. 


tions of Draco dSph would reveal the decaying dark matter 
nature of the line. The absence of the new line in stacked 
galaxy spectra of m formally excludes r DM < 1.8 x 10 28 s 
but taking into account systematical effects in spectra (e.g. 
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Projected mass density, M Sun /pc 2 

FIG. 3: Dependence of the new line flux from the expected projected 
dark matter mass. This Figure is taken from Q]; overplotted are 
the ranges for our > 2<r MOS (green) and PN (magenta) detections. 
No negative detections are shown here; the strongest restriction on 
r DM > 3.5 X 10 27 s comes from Virgo cluster, see text. 


causing significant negative residuals) and apparent uncer¬ 
tainty in dark matter distributions [17 !] produces much weaker 
bound, e.g. r DM > 3.5 x 10 27 s pD using stacked dataset 
of nearby galaxies of jl9tl with comparable exposure. Other 
bounds on decaying dark matter in ~3.55 keV energy range 
(see 119, [32,123!] and references therein) are also consistent 
with our detections after taking into account residual system¬ 
atic effects and/or uncertainties of dark matter distributions. 
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Appendix A: Dark matter distributions in galaxy clusters 


To describe dark matter distribution in galaxy clusters used in our work we compliled in Table Ullldark matter distributions from 
the literature using the extended dataset of [ 17]. 

All cluster distributions are described with Navarro-Frenk-White (NFW) profile [34]] 

p s r s 


PNFw(f) = 


r(l +r/r s ) 2 


(Al) 


parametrised by p s and r s . 


The dark matter column density inside XMM-Newton field-of-view radius R\\ = Dl x — is derived as 

2 pRia r _ 

S = —/ rdr / dz p NFW (\/r 2 + z 2 ) 

^14 JO J 


(A2) 


For the NFW density distribution <IA 1 b : 


5nfw(-R) = 


4 p s rl 

R 2 


arctan yjR 2 /r 2 — 1 


VR 2 /r 2 - 1 

Dark matter distribution parameters for our Galaxy are taken from 0. 


log 


R 

2 


(A3) 


Appendix B: Modeling astrophysical lines 


To check the astrophysical origin of the new line, we added narrow zgaussians corresponding to known astrophysical lines 
in this range. For example, according to the newest atomic database AtomDB v. 3.0.3, there are S XVI line complexes at 3.355, 
3.398, 3.424, 3.441, 3.452 and 3.460 keV. Fortunately, the intensity of these lines can be robustly predicted by the measured 
S XV line complex at ~2.63 keV, see Fig.[4]for details. 

We paid special attention to potential contribution from K XVIII lines near 3.51 keV, see e.g. [26] for details. The distance 
between these lines is smaller than the energy resolution of XMM-Newton/EPIC, so we modeled the K XVIII line complex as a 
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Object 

Reference Profile 

Rl4 

r s 

ps 

Sobj 




kpc 

kpc 

1O 6 M 0 /kpc 3 

M 0 /pc 2 

Virgo 

[35] 

NFW 

73 

560 

0.32 

808 

Centaurus 

[20 

NFW 

180 

345 

1.51 

1087 

Abell 85 

[28 

NFW 

867 

1282 

0.25 

549 

Abell 85 

[37] 

NFW 

922 

650 

0.37 

210 

Abell 478 

[38] 

NFW 

1978 

1140 

0.85 

686 

Abell 478 

□2] 

NFW 

1518 

488 

0.77 

134 

Abell 2199 

[36] 

NFW 

526 

560 

0.76 

552 

Abell 2199 

[40] 

NFW 

509 

214 

1.7 

180 

Abell 496 

[36] 

NFW 

545 

738 

0.45 

530 

Abell 496 

[37] 

NFW 

550 

420 

0.48 

191 

2A0335+096 

[28 

NFW 

593 

626 

0.52 

420 

2A0335+096 

m 

NFW 

593 

130 

3.6 

101 

Abell 1060 

m 

NFW 

195 

140 

7.2 

899 

Abell 1060 

[37] 

NFW 

211 

140 

5.8 

667 

Abell 3266 

[28 

NFW 

991 

1576 

0.19 

543 

Abell S805 

[32 

NFW 

232 

190 

2.0 

386 

Coma 

[28 

NFW 

397 

459 

1.23 

788 

Coma 

[43] 

NFW 

407 

326 

0.85 

275 

Abell S239 

[44] 

NFW 

1095 

792 

0.55 

391 

Abell S239 

[44] 

NFW 

1095 

576 

0.98 

361 

Abell 2142 

[36] 

NFW 

1469 

1654 

0.18 

406 

Abell 2142 

[23 

NFW 

1520 

990 

0.18 

144 

Abell 2319 

[36] 

NFW 

943 

1301 

0.24 

506 

Abell 1795 

[36] 

NFW 

1052 

1024 

0.34 

417 

Abell 1795 

[45] 

NFW 

1052 

393 

0.82 

139 

Abell 209 

[46] 

NFW 

3460 2513 

0.18 

408 

Abell 209 

[47J 

NFW 

3435 

502 

0.386 

24 

Perseus 

[48] 

NFW 

305 

360 

1.1 

563 

PKS 0745-191 

[28 

NFW 

1665 

1148 

0.33 

324 

PKS 0745-191 

[49J 

NFW 

1779 

230 

2.5 

59 

Triangulum 

[36] 

NFW 

856 

666 

0.83 

534 


TABLE III: Parameters of dark matter distributions of galaxy clusters used in this work. 


single zgaussian with mean energy 3.51 keV. Because there is no “reference” Potassium line to reproduce the 3.51 keV line 
flux, we fixed only the upper bound of the 3.51 keV line intensity relating it to S XVI line flux at 2.63 keV (or, if 2.63 keV is 
not detected in the dataset - to its 2er upper bound) using the procedure described in (26}]. To derive electron temperature, we 
used flux ratios of strong elemental lines, namely S XV lines at 2.45 keV, S XVI lines at 2.63 keV, Ca XIX lines at 3.90 keV and 
Ca XX lines at 4.10 keV. Because 3.51/2.63 keV line ratio is a decreasing function of electron temperature T e , see Fig. [4] we 
used minimal temperature T et u ne = min \T e< s, T e ,Ca\ for conservative estimate. Another source of uncertainty comes from the 
(largely unknown) relative K/S abundance ratio. To account possible uncertainties [ fi HH f zstl we allowed this ratio to be up to 
3 Solar values of H23ll . Note that from comparison of columns 4 and 6 one can derive that to explain new line emission solely in 
terms of K XVIII line complex at 3.51 keV, one should assume strongly supersolar (Abund[K]/Abund[S] >15 Solar) ratios for 
all our of detections. 
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FIG. 4: Line emissivity ratios for K XVIII line complex at ~3.51 keV and other emission lines of our interest as functions of the electron 
temperature T e in the plasma. The line emissivities are calculated using AtomDB version 3.0.3 with line emissivities > 10“ 22 ph cm 3 /s. 










